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ABSTRACT 

In  recent  years,  many  researchers  have  used  numerical  simulation  to  investigate  scramjet  engine  because, 
setting  up  of  testing  facilities  for  scramjet  engine  is  very  expensive.  Completion  of  simulation  work  in  less  time  is  cost- 
effective  but,  it  is  also  important  to  achieve  good  outcome  from  the  simulation  work.  In  this  context,  the  use  of 
appropriate  turbulent  model  for  the  simulation  for  DLR  scramjet  engine  is  a  vital  aspect.  From  the  available  literature, 
it  is  seen  that  many  researchers  have  used  Reynolds  stress,  SST  k-  co  and  RNG  k-s  models  in  the  field  of  supersonic  flow 
and  combustion.  Thus,  in  this  study,  the  primary  objective  was  to  find  the  better  turbulence  model  in  between  Reynolds 
stress,  SST  k-  co  and  RNG  k-e  models  for  numerical  simulation  ofDLR  scramjet  engine.  Finally,  from  the  study,  it  is  found 
that  using  SST  k-  co  model  for  numerical  simulation  of  DLR  scramjet  engine  results  in  more  accuracy  at  a  lesser 
computational  time. 
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INTRODU  CTION 

Need  of  the  modern  age  is  to  develop  reliable  and  efficient  high-speed  air  breathing  propulsion  system 
(i.e.  supersonic  to  hypersonic  flight  regime),  to  achieve  the  goal  towards  development  of  space  industries,  defence 
industries  and  easily  affordable  long-range  human  transportation  system  [1].  Turbofan  engines  are  capable  of 
operating  up  to  a  speed  of  Mach  3  and  beyond  this  Mach;  turbofan  engines  lose  its  functionality  [2].  Even,  ramjet 
propulsion  systems  are  also  not  capable  of  operating  beyond  Mach  5.  It  is  due  to  excess  pressure  and  temperature 
losses  resulting  from  deceleration  of  supersonic  flow  to  subsonic  inside  combustion  chamber  [2].  In  this  context, 
scramjet  is  one  of  the  best  available  systems  for  high  speed  propulsion  even  beyond  Mach  5 [3].  Development  of 
scramjet  engine  is  very  challenging  for  researchers  due  to  its  complex  flow  phenomena  and  unstable  chemical 
reaction  process  in  the  combustion  chamber.  Besides  this,  setting  up  of  testing  facilities  for  scramjet  engine  is  very 
expensive  [2].  The  best  economical  solution  for  investigation  of  scramjet  engine  has  attracted  the  researchers 
throughout  the  world  towards  computational  fluid  dynamics  (CFD). 

In  order  to  further  develop  the  scramjet  engine,  various  researchers  have  implemented  several  fuel 
injection  technics  [4].  One  of  the  popular  injection  schemes  for  scramjet  engine  is  transverse  injection  of  fuel  into 
supersonic  cross  flow  [5,  6].  This  results  in  formation  of  bow  shock  which  is  normal  to  the  supersonic  cross  flow. 
This  phenomenon  increases  the  pressure  loss  significantly  in  the  engine,  which  effects  its  performance.  However,  it 
is  found  that  bow  shock  wave  increases  combustion  performance  by  creating  a  separation  zone  at  downstream  of 
the  shock  which  supports  the  holding  of  the  flame  [7].  Researchers  in  the  recent  past  have  suggested  the  use  of 
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parallel  fuel  injection  in  scramjet  engine  by  using  strut-based  injector  to  inject  fuel  into  the  core  of  the  supersonic  fiow. 
This  results  in  reduction  of  the  effect  of  normal  bow  shock  that  minimize  the  performance  of  engine  [8,  9].  However,  strut- 
based  injection  also  results  in  slight  drop  in  mixing  performance  as  compared  with  transverse  injection  [8,  9]. 

Development  of  various  numerical  modelling  in  computational  fluid  dynamics  plays  significant  role  in 
investigating  supersonic  combustion  precisely  with  considerable  computational  cost.  Conventionally,  developed  Reynolds 
Averaged  Navier  Stokes  (RANS)  models  are  not  enough  to  analyse  supersonic  combustion  accurately.  LES  (Large  eddy 
simulation)  model  has  a  significant  achievement  in  the  field  of  numerical  modelling  of  supersonic  combustion,  which  is 
capable  of  predicting  quantitative  and  qualitative  outcome  of  scramjet  engine  [10,  11].  LES  model  has  been  introduced  by 
many  scientists,  in  the  process  of  investigating  the  operation  of  strut  scramjet  engine.  It  is  seen  that  using  LES  model  for 
numerical  simulation  results  in  reasonably  high  computational  cost  and  time  [1,  2].  However,  in  recent  years,  many 
researchers  have  used  other  computational  models  to  investigate  scramjet  engine.  Dinde  et  al.  [12]  have  numerically 
investigated  scramjet  combustor  by  introducing  RNG  k -s  turbulence  model  with  finite  rate  single  step  chemistry  model. 
The  researchers  found  convincing  results  of  scramjet  combustor  in  regard  to  H2  combustion.  Later,  Luo  et  al.  [13]  conducted  a 
numerical  investigation  on  struct  scramjet  combustor  and  found  that  SST  k-  co  model  is  less  convincing  compared  with  k -s  for 
non-reacting  flow  simulation.  SST  k-  co  model  introduced  by  Pandey  et  al.  [14]  in  their  simulation  work  of  scramjet  combustor 
observed  that  performance  of  scramjet  engine  is  dependent  on  effective  mixing  of  fuel  into  air  and  reduction  of  pressure  loses. 
The  result  obtained  by  Pandey  et  al.  using  SST  k-  co  model  were  highly  accurate  in  comparison  with  available  experimental 
data.  Huang  et  al.  [15]  numerically  examined  effects  of  parametric  changes  on  strut  scramjet  combustion  flow  field,  and  it  is 
highlighted  in  their  research  that  two  step  reaction  mechanism  with  SST  k-  co  model  has  slight  changes  with  single  step 
mechanism.  From  their  investigation,  it  also  can  be  seen  that  on  using  SST  k-  co  model,  the  results  obtained  are  very  close  to  the 
available  experimental  values.  Reynolds  stress  model  for  supersonic  flow  simulation  was  conducted  by  Frauholz  et  al.  [16],  and 
they  concluded  that  the  model  is  capable  of  predicting  boundary  layer  shockwave  interaction  effectively  and  effect  of  grid 
reduction  is  invariant  of  solution  outcome.  However,  in  the  numerical  simulation  of  strut  scramjet  engine,  Reynolds  Stress 
model  has  not  yet  been  used.  Hence,  the  primary  objective  of  this  paper  is  to  find  the  better  turbulence  model  in  between 
Reynolds  stress,  SST  k-  co  and  RNG  k -s  models  for  numerical  simulation  of  DLR  scramjet  engine.  Moreover,  in  this  paper,  the 
results  of  all  the  numerical  simulation  is  compare  against  the  experimental  result  obtained  by  Waidmann  et  al.  [17-20]  of 
German  aerospace  centre  (DLR). 

COMPUTATIONAL  DOMAIN  AND  GRID  GENERATION 

Completion  of  simulation  work  in  less  time  is  cost-effective  but,  it  is  also  important  to  achieve  good  outcome 
from  the  simulation  work.  So,  in  our  present  work  we  have  used  2D  models  with  structured  mesh  for  simulation.  To  design 
the  2D  model  and  to  generate  the  structured  mesh,  ICEM  of  Ansys  software  is  used.  In  this  work,  the  dimensions  of  2D 
model  are  kept  similar  to  the  experimental  model  by  waidmann  et  al.  [17-20]  as  shown  in  Figure  l(a).  The  2D  model  has  a 
length  of  340  mm  with  upper  wall  diverging  at  3°  in  outward  direction.  The  distance  between  upper  and  lower  wall  at  inlet 
is  50  mm.  At  the  centre  of  the  combustion  chamber,  a  strut  of  32mm  length  and  12°  strut  angel  is  positioned  at  77mm  from 
the  start  of  the  Scramjet  engine.  Moreover,  structured  grid  size  of  325200  elements  and  Y+  value  of  less  than  2.5  near  the 
walls  is  considered  for  simulation  work,  which  is  shown  in  Figure  l(b)  [21].  For  grid  dependency  analysis,  results  of 
bottom  wall  pressure  within  Scramjet  combustor  for  non-reacting  flow  are  observed  which  is  shown  in  Figure  2.  Primarily, 
three  different  grid  sizes  of  122980,  226340  and  325200  are  used  in  this  analysis  for  three  models.  In  case  of  Reynolds 
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Stress  model  analysis,  one  more  grid  size  of  424412  is  added,  as  it  is  observed  from  Figure  2(c)  that  there  is  a  significant 
deviation  between  obtained  results  from  grid  size  226340  and  325200.  Comparing  grid  size  of  325200  with  others  for  all 
three  models,  no  noteworthy  change  in  the  simulation  results  are  obtained(i.e.  grid  independent  solution  is  obtained)[22]. 
In  the  generated  grid  inside  the  model,  all  walls  are  densely  meshed  in  order  to  properly  calculate  shockwave  boundary 
layer  interaction.  Further,  dense  grid  size  is  also  maintained  at  inlets,  outlet  and  shear  layers.  In  Figure  l(b),  the  grid  pattern 
that  is  shown  has  half  of  the  original  grid  size  for  better  realization. 


(a)  (b) 

Figure  1:  (a)  2D  Diagram  of  DLR  Scramjet  Engine  (b)  Grid  Pattern  for  2D  DLR  Model 
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Figure  2:  Grid  Dependency  Test  using  Bottom  Wall  Pressure  for 
(a)  RNG  K-£,  (b)  SST  K-O  and  (c)  Reynolds  Stress  Model 

NUMERICAL  MODELLING  AND  GOVERNING  EQUATIONS 


In  this  study,  all  the  2D  simulations  are  conducted  at  steady  state  condition  in  Ansys  Fluent  software.  RNG  k-e , 
SST  k-co  and  Reynolds  Stress  model  are  the  three  turbulence  models,  which  are  selected  for  the  present  numerical 
investigation.  For  combustion  modelling,  hybrid  chemistry  model  of  finite  rate  and  eddy  dissipation  are  introduced  along 
with  the  mentioned  turbulence  models.  Single  step  chemical  reaction  between  hydrogen  and  oxygen  (2H2  +  02  =  2H20)  is 
considered  in  the  present  study,  to  reduce  the  computational  time  which  is  reported  by  Choubey  and  Pandey  [21]. 
Description  of  Transport  Equation  for  all  turbulence  model  used  in  this  work  is  given  below. 
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RNG  k-e  Model 

d/8t(pk)  +  d/dXiipkUi)  =  dldxj  (akpeff  +  Gk  +  Gb  -  pe  -  YM  +  Sk  (1) 

And, 

Tt  (P£)  +  £-  (P£UJ  =  £:  (<W/  +  Cu  I (Gfc  +  -  C^P  7-Rc  +  se  (2) 

Where, 

Gk,  Gb  is  the  turbulence  kinetic  energy  produced  due  to  mean  velocity  gradient  and  buoyancy,  respectively.  Ym  is 
the  ratio  of  fluctuating  dilation  to  the  overall  dissipation  rate  in  compressible  turbulence.  The  two  terms,  a£  and ak  are 
inverse  effective  Prandtl  numbers.  Further,  the  user  defined  terms  are  represented  as  S£,  Sk. 

SST  k-(o  Model 

h  0>‘) + 1,  =  ir,  (r»  £;)  +  G‘  -  r‘ + ® 

and, 

k m 4  &>“”/)  =  sj  (r»  w)  + -  Y« +  D-  * s-  +> 

Where,  G^  is  the  representation  of  oo.  Tk  and  is  the  effective  diffusivity  of  k  and  oo  respectively.  Dissipation  of 
k,od  is  represented  by  Yk,  Y^.  Further,  user  defined  terms  are  Sk  and  S^. 

Reynolds  Stress  Model 


d  f  _ ..  d  f  .  d  r  -  - 7 - rn  d  r  d  f _ / - dU:  - dUj\ 

+  =  -7Tklpu'^  +  p,(^u;  +  5-^')J +  a^fe(u^)J  -  du;u^  +  u;u^J- 

pPidiufi  +  9j<e)  +  p'  -  2pnk(UJu£eikm  +  <ujjejkm)  + 


(5) 


BOUNDARY  CONDITIONS 

In  this  study,  supersonic  conditions  are  applied  at  both  air  and  hydrogen  inlet,  which  is  incorporated  by  pressure 
far-field  condition  in  Ansys  fiuent  software.  For  all  the  simulation,  pressure  outlet  boundary  condition  is  applied  at  the 
outlet  of  the  model[23].  All  fixed  wall  including  strut  is  introduced  with  no  slip  condition  for  the  simulations. 
All  the  values  of  variable  at  inlet  for  air  and  hydrogen  are  given  in  Table  1.  The  convergence  criteria  are  considered  on  the 
basic  of  mass  imbalance  test,  which  is  the  difference  between  inflow  mass  and  out  flow  mass  of  computational  domain 
(min  -  niout)  [24].  Hence,  for  all  the  simulations,  mass  imbalance  up  to  0.001  is  set  as  a  criterion  for  convergence. 


Table  1:  Inflow  Parameters  for  Numerical 
Simulation  of  DLR  Scramjet  Engine 


Variables 

U  (m/s) 

Ma 

T(K) 

P(Pa) 

Yo2 

Yr2 

Yh20 

YN2 

Air 

730 

2 

340 

100000 

0.232 

0 

0.032 

0.736 

Hydrogen 

1200 

1 

232 

100000 

0 

0 

0 

1 
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Experimental  results,  observed  by  Waidmann  et  al.  [17-20]  are  taken  as  reference  to  validate  the  present  simulation 
work.  Figure  3(a)  shows  the  Schlieren  image  of  DLR  experiment  and  Figures  3(b),  3(c)  and  3(d)display  the  Mach  no  contour  of 
cold  flow  simulation  using  RNG  k-e,  SST  k-co  and  Reynolds  Stress  model,  respectively. 


Figure  3:  (a)  Schlieren  Image  (b)  RNG  K-£  Model  Predicted  Mach  No  Contour  (c)SST  K-O 
Model  Predicted  Mach  No  Contour  (d)  Reynolds  Stress  Model  Predicted  Mach  No  Contour 


(a)  (b)  (c) 

Figure  4:  (a)  Distribution  of  Bottom  Wall  Pressure  of  DLR  Experiment  [34-36] 
with  Numerical  Calculation  (b)  Plot  of  Percentage  Deviation 
for  Bottom  Wall  Pressure  of  three  Models 


From  figures,  it  is  clearly  seen  that  the  shock  structure  predicted  by  the  simulations  and  Schlieren  image  of  DLR 
are  similar  in  structure.  Moreover,  a  comparison  between  bottom  wall  pressure  distribution  obtained  from  the  literature  of 
Waidmann  et  al.  [17-20]  and  current  simulations  is  also  made  in  this  paper  which  is  highlighted  in  Figure  4.  From  Figure 
4(a),  it  is  observed  that  pressure  distribution  predicted  along  the  bottom  wall  by  SST k-co  model  is  nearest  to  the 
experimental  results  available  in  the  literature  [17-20]  in  comparison  with  other  turbulence  models  used  in  the  study.  In 
addition  to  the  above,  Figure  4(b)  shows  the  percentage  change  in  prediction  of  pressure  along  the  bottom  wall  of  DLR 
scramjet  engine  by  SST k-co,  RNG  k-s  and  Reynolds  Stress  models  in  comparison  with  the  experimental  results  available  in 
the  literature[  17-20]. 

RESULTS  AND  DISCUSSIONS 


In  this  study,  both  cold  flows  as  well  as  reactive  flow  simulations  are  conducted  using  RNG  k-s,  SST  k-co  and 
Reynolds  Stress  model.  In  case  of  simulations  of  reactive  flow,  an  equivalence  ratio  of  0.136is  maintained  which  is  equal 
to  the  equivalence  ratio  used  by  Waidmann  et  al.  [17-20]  in  their  experimental  investigation  of  DLR  scramjet  engine.  Solution 
scheme  used  in  all  the  numerical  simulations  is  second  order  upwind,  which  is  highly  efficient  for  capturing  shockwave 
structure  and  boundary  layer  accurately. 
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Contour  of  Mach  No.  and  static  temperature  for  different  turbulence  models  is  shown  in  Figure  5.  On  comparing 
Mach  no  and  static  temperature  contour  for  RNG  k-s  and  SST  k-co  models,  it  is  observed  that  there  is  no  significant 
difference.  Both  the  models  predict  same  static  temperature.  However,  on  using  Reynolds  Stress  model,  shock  structure 
(in  Mach  no  contour)  and  temperature  contour  of  different  nature  are  found  as  compare  to  RNG  k-s  and  SST  k-co  models. 
Simulation  using  Reynolds  Stress  model  has  predicted  highest  temperature  of  2353K,  which  is  approximately  14%  higher 
than  the  highest  temperature  predicted  by  RNG  k-s  and  SST  k-co  models.  Reynolds  stress  model  estimates  high  static 
temperature,  because  it  predicts  the  formation  of  eddies  in  few  regions  at  the  downstream  of  the  flow  field  which  act  as 
flame  stabilizer. 


Figure  5:  Mach  No  and  Static  Temperature  Contour  For  RNG  K-E  (a,  b), 
SST  K-O  (c,  d)  and  Reynolds  Stress(e,  f)  Models 


Mixing  efficiency  is  a  very  important  parameter  for  the  designing  of  scramjet  engine  combustor.  Good  mixing  of 
air  and  fuel  in  a  short  combustor  distance  provides  shorter  combustor  length,  which  results  in  reduction  of  weight  of  the 
combustor  [21].  Mixing  efficiency  at  section  X,  for  a  non-reacting  flow  is  given  by  [25]. 


Vmix(X) 


f  apuYH2  dA 

mH2  00 


(6) 


Where,  Area  of  channel  cross  section  is  A,  mH2(X )  is  mass  flow  of  hydrogen  at  section  X  and a  is  given  by 


a 


Here,  6  is  equivalence  ratio. 


Figure  6(a),  shows  mixing  efficiency  plot  along  the  combustor  chamber  for  different  models.  Out  of  these  models, 
RNG  k-s  takes  shortest  distance(163mm)  for  complete  air-fuel  mixing  and  Reynolds  Stress  model  takes  longer  distance 
compared  with  all  other  models.  Moreover,  Figure  6(b)  shows  the  percentage  change  in  prediction  of  mixing  efficiency  of 
DLR  scramjet  engine  by  RNG  k-s  and  Reynolds  Stress  models  in  comparison  with  SST  k-co  model.  The  reason  for 
comparing  the  results  of  RNG  k-s  and  Reynolds  Stress  models  with  that  of  the  SST k-co  model  is  because,  the  results 
predicted  by  SST k-co  model  is  found  to  be  more  accurate  in  term  of  experimental  results  which  are  available  in  the  existing 
literature.  From  Figure  8,  it  is  seen  that  the  average  deviation  in  mixing  efficiency  predicted  by  RNG  k-s  model  is  less  than 
Reynolds  Stress  model.  For  Reynolds  Stress  model,  higher  deviation  of  13.85%  is  observed  at  a  combustor  length  of 
166mm  and  that  for  RNG  k-s  model ,  higher  deviation  of  23.27%  is  observed  at  a  combustor  length  of  157mm. 
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(a)  (b) 

Figure  6:  (a)  Mixing  Efficiency  (r|mix)  Plot  For  RNG  K-E,  SST  K-O  And  Reynolds  Stress  Models 
(b)  %  Deviation  In  Mixing  Efficiency  For  RNG  K-£  And  Reynolds  Stress 
Models  In  Comparison  With  SST  K-O  Model 

Another  important  parameter  in  a  scramjet  engine  is  combustion  efficiency.  Study  of  combustion  phenomena  in  a 
scramjet  engine  is  very  important  to  finding  out  the  strength  of  the  scramjet  engine,  and  it  is  determined  by  combustion 
efficiency  (?7C)[21]. 


1  _  /  A(x)pgasuYH2  dA 
mH2  (inj) 


(7) 


where,  mass  flow  rate  of  injected  hydrogen  is  mH2^  y 


Figure  7(a),  shows  combustion  efficiency  plot  along  the  combustor  chamber  for  different  models.  Out  of  these 
models,  RNG  k-s  and  SST  k-co  predicted  similar  nature  of  combustion  efficiency  up  to  combustor  length  of  140  mm. 
Beyond  the  140mm  length  of  combustor,  RNG  k-e  model  predicted  higher  combustion  efficiency  than  SST  k-co  model. 
RNG  k-s  takes  shortest  distance(163mm)  for  complete  air-fuel  mixing  and  Reynolds  Stress  model  takes  longer  distance 
compared  with  all  other  models  in  this  study.  Furthermore,  Figure  7(b)  shows  the  percentage  change  in  prediction  of 
combustion  efficiency  of  DLR  scramjet  engine  by  RNG  k-s  and  Reynolds  Stress  models  in  comparison  with 
SST  k -co  model. 
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Figure  7:  (a)  Combustion  Performance  Plot  between  RNG  k-r,  SST  k-co  and  Reynolds  Stress 
Models  (b)%  Deviation  in  Combustion  Efficiency  for  RNG  k-r  and  Reynolds 
Stress  Models  in  Comparison  with  SST  k-co  Model 

The  purpose  for  comparing  the  results  of  RNG  k-s  and  Reynolds  Stress  models  with  that  of  the  SST k-co  model  is  for 
the  reason  that  the  results  projected  by  SST k-co  model  is  found  to  be  more  precise  in  terms  of  experimental  results,  which 
are  presented  in  the  available  literature.  From  Figure  10,  it  is  seen  that  the  average  deviation  in  combustion  efficiency 
predicted  by  RNG  k-s  model  is  less  compared  with  Reynolds  Stress  model.  For  Reynolds  Stress  model,  higher  deviation  of 
26.05%  is  observed  at  a  combustor  length  of  210mm  and  that  for  RNG  k-s  model ,  higher  deviation  of  5.65%  is  observed  at 
a  combustor  length  of  166mm. 
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The  simulation  in  this  study  is  conducted  on  a  workstation  equipped  with  Intel  Xenon  processor  of  8  GB  memory. 
The  computational  time  spend  for  cold  flow  simulation  using  RNG  k-s  and  SST  k-co  models  were  approximately  16  hours 
and  that  for  using  Reynolds  Stress  model  was  31  hours.  Whereas,  the  computational  time  spends  for  reactive  flow 
simulation  using  RNG  k-s  and  SST  k-co  models  were  approximately  36.5  hours  and  that  for  using  Reynolds  Stress  model 
was  approximately  74  hours. 

CONCLUSIONS 

In  this  study,  the  primary  objective  was  to  find  the  better  turbulence  model  in  between  Reynolds  stress,  SST  k-  co  and 
RNG  k -s  models  for  numerical  simulation  of  DLR  scramjet  engine.  Base  on  the  results  found  from  the  simulations,  the 
following  conclusions  are  drawn. 

•  The  results  predicted  by  SST k-co  model  is  found  to  be  more  accurate  than  the  results  predicted  by  RNG  k-s  and 
Reynolds  Stress  models  in  comparison  with  experimental  results. 

•  In  reference  to  the  results  obtained  by  using  SST k-co  model ,  it  is  seen  that  the  average  deviation  in  mixing  and 
combustion  efficiency  predicted  by  RNG  k-s  model  is  less  compared  with  Reynolds  Stress  model. 

•  The  computational  time  required  for  simulations  using  RNG  k-s  and  SST  k-co  models  are  substantially  less  in 
comparison  with  Reynolds  Stress  model. 

•  Thus,  from  this  study,  it  can  be  finally  establish  that  the  use  of  SST  k-  co  model  for  numerical  simulation  of  DLR 
scramjet  engine  will  provide  more  accurate  result  at  a  lesser  computational  time. 
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TERMINOLOGY 

•  DLR  -  German  Aerospace  Centre  (Deutsches  Zentrum  fur  Luft-und  Raumfahrte.  V.). 

•  2D  -  Two-dimensional,  k  -  Turbulent  kinetic  energy,  SST  -  Shear  stress  transport,  co  -  Turbulent  dissipation  rate. 
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